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Abstract. Univariate spline discrete quasi-interpolants (abbr. dQIs) 
are approximation operators using B-spline expansions with coeffi- 
cients which are linear combinations of discrete values of the function 
to be approximated. When working with nonuniform partitions, the 
main challenge is to find dQIs which have both good approximation 
orders and bounded uniform norms independent of the given partition. 
Near-best dQIs are obtained by minimizing an upper bound of the in- 
finite norm of dQIs depending on a certain number of free parameters, 
thus reducing this norm. This paper is devoted to the study of some 
families of near-best dQIs of approximation order 2. 



§1. Introduction 

A spline quasi-interpolant (abbr. QI) of / has the general form 

Q/= $> Q (/)B a 

aeA 

where {B a , a G A} is a family of B-splines forming a partition of unity and 
(/)) a ^} i s a family of linear functionals which are local in the sense 
that they only use values of / in some neighbourhood of E a = supp(B a ). 
The main interest of QIs is that they provide good approximants of func- 
tions without solving any linear system of equations. In the literature, one 
can find the three following types of QIs: 

(i) Differential QIs (abbr. DQIs) : the linear functionals are linear combi- 
nations of values of derivatives of / at some point in E a (see e.g. [5-7]). 

(ii) Discrete QIs (abbr. dQIs) : the linear functionals are linear combina- 
tions of values of / at some points in some neighbourhood of E Q (see e.g. 
[1-3], [6], [9], [11], [13], [15-16], [24]). 

(iii) Integral QIs (abbr. iQIs) : the linear functionals are linear combina- 
tions of weighted mean values of / in some neighbourhood of S a (see e.g. 
[2-3], [6], [13-14], [24-25]). 
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In this paper and a subsequent one, we shall study various types of uni- 
varaite dQIs and iQIs, more specifically those that we call near-best QIs 
which are defined as follows: 

(dQIs) assume that /U Q (/) = J2 /3eFa ^ a (P)f( x p) where the finite set of 
points {xp,(3 G F a } lies in some neighbourhood of S a . Then it is clear 
that, for H/lloo < 1 an d a G A, |/i a (/)| < ||A a ||i, where A a is the vector 
with components A a (/3), from which we deduce immediately 

IIQIloo < Y] \fi a (f)\B a < max|// a (/)| < max || A a ||i = ui(Q). 

aeA 

Now, assuming that n = card(F a ) for all a, we can try to find a A* G M n 
solution of the minimization problem 

||A* ||i = min{||A Q ||i; A Q G IT, V a X a = b a } 

where the linear constraints express that Q is exact on some subspace of 
polynomials. Thus, we finally obtain 

IIQIloo < v{{Q) = max || A a ||i. 

aeA 

(iQIs) assume that /J, a (f) = J2(3eF a ^a(P) / s Mp(t)f(t)dt, where the B- 
splines Mp are normalized by J Mp = 1. Once again, for ||/||oo < 1) w ^ 
have 

whence, as we obtained above for dQIs, 

IIQIloo < max ||A Q ||i = v{(Q). 

ct£A 

As emphasized by de Boor (see e.g. [5], chapter XII), a QI defined on 
non uniform partitions has to be uniformly bounded independently of the 
partition (abbr. UB) in order to be interesting for applications. Therefore, 
the aim of this paper is to define some families of dQIs satisfying this 
property and having the smallest possible norm. As in general it is difficult 
to minimize the true norm of the operator, we have chosen to solve the 
minimization problems defined above. A further paper [25] will develop 
the case of iQIs on nonuniform partitions. A few results are given in [24]. 

The paper extends some results of [1] [13] and is organized as follows. We 
first recall some "classical" QIs of various types and we verify that they 
are UB. Then we define and study several families of discrete and integral 
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QIs, depending on a finite number of parameters, for which we can find 
i>i(Q). We show that this problem has always a solution (in general non 
unique). We give more specific examples for quadratic and cubic splines. 
Of particular interest are the results of theorems 3,5 and 6 where we show 
that some families of dQIs are uniformly bounded independently of the 
partition. Finally, we briefly give some applications to the approximation 
of functions, to quadrature formulas and to pseudo-spectral methods (see 
e.g. [12], [29]). A parallel study of spline QIs is done in [2] for uniform 
partitions of the real line and in [3] for some uniform triangulations of the 
plane. 



§2. Notations 

We shall use classical B-splines of degree m on a bounded interval / = [a, b] 
or on / = M. For the sake of simplicity, in the case / = M, we take an 
increasing sequence of knots T = G 7Z}. In the case / = [a, 6], we 
take the usual sequence T of knots defined by (see [5] [11] [19] [28]): 

a, = t— m = . . . = to, b = t n = . . . = t n _|_ m 

a < ti < t 2 < ■ ■ ■ < t n -i < b 

For J = {0, . . . , n + m — 1}, the family of B-splines {Bj,j G J}, with 
support Ej = [tj- m ,tj + i] is a basis of the space S m (I,T) of splines of 
degree m on the interval / endowed with the partition T. These B-splines 
form a partition of unity, i.e. J^jeJ^j = 1- We denote hi = ti — £j_i for 
all indices i. 

Let JM m =Wn[0,m-l] and Tj = {tj_ r ,r G N m }: we recall that the 
elementary symmetric functions cr;(T) of the m variables in Tj are defined 
by <Jo(Tj) = and for 1 < I < m, by 

a l(Tj) = tj- ri tj- r2 ...tj- rr 
0<ri<r2<---<ri<m — l 

Let C l m = be the binomial coefficients, then the monomials ei (x) = 

x l can be written e { = J2 ie j d i ^ with ^ = a i( T t)/ c L for < / < m. 
This is a direct consequence of Marsden's identity ([3], chapter IX). 

§3. Differential QIs 

For all j G J, we define ij)j(t) = Ylrew (tj-r ~ t) (thus ifjj G P m for all 
j G J). From de Boor and Fix [7] or de Boor ([5], chapter IX), we know 
that for any r G £-,-, the functionals 

= -yE(- 1 ) m "^ m "^ J (T)^7(r) 

lib' 

1=0 
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are dual functionals of B-splines, i.e. they satisfy, for all pairs G JxJ 

Xj(Bi) = Sij. 

Therefore the differential quasi-interpolant (abbr. DQI) 

jeJ 

satisfies QBj = Bj for all j G J, i.e. Q is a projector on the space 
S m (I, T). In practice, it is interesting to choose r = 9j = ^ Z) s gK m ^i- s = 
si(Tj) = 6y. However, the computation of Xj(f) needs the evaluation all 
derivatives of polynomials ipj. Another method consists in writing Q in 
the form 

ieJ 

whose coefficient functionals are defined by 



1=0 

and to impose that Q be exact on monomials of degree at most m 

Qek = efc for < k < m. 

Setting ai(9i) = 9~ l ai(9i) and (3 s (9i) = 9~ s 9\ s \ we obtain the following 
system of linear equations, for < s < m: 



Y J C l s a l {9 i ) = (3 8 {9 i ). 



1=0 

The solutions of this system are given by 



a s {9 i ) = Y J {-^) s ~ l C%{9 i ). 



1=0 

Thus we finally obtain 

Theorem 1. The coefficients of the differential forms {Aj(/), i G J}, are 
given, for < s < m, by 



..((?<) = X)(-i)- i ciO? ) - 



1=0 
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However, these DQIs need the values of derivatives of /, so they are not 
very easy to use in applications and we will not study them any more. Let 
us only give examples of quadratic and cubic DQIs. 

Example 1: Quadratic spline DQIs (see also section 4 below). In de 
Boor's form, we have for r = 9f Xj(f) = f(9j) - \{9) - ef ) )D 2 f(6 j ), 
where 9j = \(tj-\ + tj) and 9j = tj-itj, whence 9 2 — 9^ = \{tj-\ — 
tjf and finally Xj(f) = f{9 3 ) - \h 2 D 2 f(9j). Theorem 1 gives a Q {9 j ) = 
1, ai (9j) = 0, a 2 (%) = 9) - 9f\ hence Xj{f) = Xj(f). 

Example 2: Cubic spline DQIs. In de Boor's form, we have for r = Qy. 
W) = f(9 )-\(9]-9f ) )D*f(9 J )-y 3 (9 J )D*f(9 3 ), with 0, = §(^_ 2 + 
tj-i + tj) and 6>j 2) = ±(i,_ 2 t,-i + tj-itj + tj- 2 tj), whence 9 2 - 9y = 
Kh^ + hi-ihi + h?) andVj(^) = (t j -9 J )(t J _ 1 -9 J )(t J -2-9 j ). Finally, 
we obtain Xj(f) = f(9 J ) - ^(/i 2 _! + K_ X K + h 2 )D 2 f{9 J ) - I fe(2/i i _ 1 + 
hi){hi — hi-i)(hi-i + 2hi)D^f(9j). Theorem 1 gives a (9j) = l,ai(9j) = 

0, a 2 (9j) = 9) - 9 { f\a z {9 3 ) = 0< 3) - W 9^ + 29) = ±{2K-i + K)^ - 
hi-x^hi-x + 2hi), whence Xj(f) = Xj(f). 

§4. Uniformly bounded discrete QIs exact on P 2 

It is now possible to derive discrete QIs from the preceding DQIs by re- 
placing the values of derivatives D f{6j)/l\ of / by divided differences at 
the points 9 r lying in Sj. Doing this, we loose the property of projection 
on S m (I,T). However, by choosing conveniently the divided differences, 
we can obtain some families of dQIs which are UB and exact on specific 
subspaces of polynomials. 

Let us construct for example a family of dQIs of degree m which are 
exact on P 2 . We start from functionals which are truncations at order 2 
of those of the preceding section: 

A f(/) = ^^- l ) m ~ lDm ~ l ^) Dl f^ 

lib' 

1=0 

As ipj(t) is of degree m, we obtain successively D m ipj(T) = (— l) m m!, 
D m_1 'i/>j(T) = (-l) m ml(r-9 J ), D m - 2 ^-(r) = ±(-l) m ml(T 2 -29 J t+9 j 2) ). 
More specifically, taking r = 9j, we get 

D m -Vi(*i) = 0, D™~ 2 ^{9 3 ) = ^(-irm!(#f - 9 2 ) 
and we obtain the DQI exact on P 2 

jeJ 
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whose coefficient functionals are given by 

><f\f) = f{e j )-\{e 2 j -ef ) )D*f{e j ). 

We recall the expansion (se e.g. [9] [17]): 

df^ = 6 2 — 6f^ = — — (tj-r - tj-s) 2 > 0. 

(r, S )6N^,r< S 

On the other hand, \D 2 f(9j) coincide on the space P 2 with the second 
order divided difference [9j-i,6j,6j+i]f, therefore the dQI defined by 

jeJ 

with coefficient functionals 

is also exact on P 2 . Moreover, one can write 

vf\f) = a-ifi-i + kfo + Cif i+1 
with a z = -^ (2) /A^-i(A^_i + A0j), b % = 1 + #f VA^-iA^, 

—(2) 

Ci = —6\ / A6i(A6i-i + A6{). So, according to the introduction 

\\Q* 2 \\oo <max(\ai\ + \bi\ + \ci\) < 1 + 2 max^ 2 } /A^-iA^. 

The following theorem extends a result given for quadratic splines in 
[13] [22] [23]. 

Theorem 2. For any degree m, the dQIs Q\ are uniformly bounded. 
More specifically, for all partitions of I: 

Halloo <[^(m + 4)] 

Proof: We only give the proof for m = 2k + 1, the case m = 2k being 
similar. For the sake of simplicity, we take j = k, i.e. we shall determine 
an upper bound of the ratio 

N k /D k = e% ) /Ad k _ 1 Ad k 
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with 

m z (m — 1) , , 



<r<s<m 



Setting H = Y^Li hi, then we get a lower bound for the denominator 
1 1 if 2 

D k = -J (* m -l - *0)(*m " *l) = -o (^0 + H)(H + h^) > 

m z m z m z 

The numerator iVfc is composed of k — 1 pairs of sums (S P7 S' p ) 

s p = (t r — t s ) , ^ = (^r _ ^s) 5 

s—r=p s — r —k+p—l 

for 1 < p < fe — 1. Both sums contain at most p times the terms h 2 and 
2/ij/ij (i^ j), hence we can write S p + S' < 2pH 2 , which implies 

2H 2 kS 2 

N k < -(l + 2 + ...fc-l) = 

(m — l)-%m — 2) 2(m — l) z 

so, we get 

iVfc/Ac < k/2, 
and finally, for m = 2k + 1 odd 

HQ2II00 < k + 2 = l(m + 3) = [i(m + 4)] 

For m = 2/c, we obtain respectively D k > jg^ and Nk < 4 ^ 2k -i) ■> whence 
Nk/Dk < ^rzj, an d finally for m = 2k even 

\\Q* 2 \\<k + 2 = ^(m + 4) = {±(m + 4)] 

□ 



§6. Existence and characterization of near-best discrete QIs 
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6.1. Existence of near-best dQIs 

We consider the following family of dQIs defined, for the sake of simplicity, 
on / = 1R endowed with an arbitrary non-uniform increasing sequence of 
knots T = {U;i E 7L}, 

Qf = Q p , q f = J2^(f) B i- 
Their coefficient functionals depend on 2p + 1 parameters, with p > m 

s =—p 

and they are exact on the space P g , where q < min(m, 2p). The latter 
condition is equivalent to Qe r = e r for all monomials of degrees < r < q. 
It implies that for all indices i, the parameters Xi(s) satisfy the system of 
q + 1 linear equations: 

X i(s)Q r i+ s = 0i r \ 0<r<q. 

s=—p 

The matrix V { E ^ q+1)x{2p+1) of this system, with coefficients V^r, s) = 
0l +s , is a Vandermonde matrix of maximal rank q + 1, therefore there are 
2p — q free parameters. Denoting hi E 1R <?+1 the vector in the right hand 
side, with components bi{r) = 6*^, < r < q, we consider the sequence 
of minimization problems, for i E TL\ 

min || Aj||i, Vi\ = hi. 

We have seen in the introduction that z^i(<5) = maxi 6 x; min ||Ai||i is an 
upper bound of ||<2q||oo which is easier to evaluate than the true norm of 
the dQI. 

Theorem 3. The above minimization problems have always solutions, 
which, in general, are non unique. 

Proof: The objective function being convex and the domains being affine 
subspaces, these classical optimization problems have always solutions, in 
general non unique. 

We postpone to sections 7 and 8 the computation of some optimal solutions 
in the case q = 2. 
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6.2. Characterization of optimal solutions 

For b G R m and A e R mXn , let us consider the /i-minimization problem 

(1) min ||r(a)||i, r(a) = b — Aa. 

We recall the characterization of optimal solutions for Zi-problems given 
in [30], chapter 6. Define the sets 

Z(a) = {1 < i < m\ r»(a) = 0} 

V(a) = {ve R m ; ||v||oo < 1, «* = sgn{r l {a)) for i <£ Z(a)} 

Theorem 4. a* is a solution of (1 ) if and only if there exists a vector 
v* e V(a*) satisfying A T v* = 0. 

§7. A general family of spline discrete QIs exact on P 2 

In this section, we restrict our study to the subfamily of spline dQIs which 
are exact on P 2 , i.e; we choose q = 2. We shall try to characterize optimal 
solutions in the sense of theorem 3 with the help of theorem 4. Let 

Q P ,2f = 

where the coefficient functionals depend on 2p + 1 parameters 

M = E Mr)/(0 i+r ). 

r=—p 

We shall need the following sets of indices 

K = {-p,...,p}, iT = {-p,0,p}, K = K\K*. 

K = K 1 UK 2 , K 1 = {-p+l,...,-l}, K 2 = {l,...,p-1}. 
The three equations expressing the exactness of Q Pj2 on P 2 can be written 

Ai(-p) + Ai(0) + Kip) = 1 - A '( r ) 

reK 

Bt-pXii-p) + 9M0) + e i+p \i(p) = 9i - E e r \i{r) 

reK 
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This system has a positive Vandermonde determinant 

Vi = V($i-p, 9i, i+p ) = ($i — 9i- p )(6i +p — 6i)(9 i+p — 9i-p). 
Let us denote by 

(A*(-p),A*(0)A*(p)) 

the unique solution of the above system with the right-hand side obtained 
by taking Aj(r) = for all r £ K. Using Cramer's rule and the determi- 
nants Wi(s) obtained by replacing the column of 9i +s in Vj by this rhs, 
we obtain 

K(-p) = Wi(- P )/Vi, a*(o) = Wi(o)/Vi, \*(p) = Wi(p)/Vi. 

Then we can express the general solution of the above system in the form 

Ai(0) = A*(0) - Pr\{r) - PMs) 

\{P) = Kip) + Yl 7rA l (r) - Y ls\(s) 
The various coefficients are quotients of Vandermonde determinants 

a r = v(e r , 0i, e i+p )/Vi, a s = v{6i, e s , e i+p )/Vi, 

Pr = V(9i- P , 9 r , 9 i+p )/Vi, p s = V(9i- P , 9 S , 9 i+p )/Vi, 
lr = V(6i- P , 9 r , 6i)/Vi, 7s = V(9 i _ p , 9,, 9 S )/V i . 
We denote by Q* 2 the spline dQI whose coefficient functionals are 

(/) = Ki-P)f{0i- P ) + A?(0)/(tfi) + X*(p)f(9 i+P ). 

In that case, an upper bound of the norm is max i€ ^ v* where 

^*HA*(-p)| + |A*(0)| + |A*(p)| 

Theorem 5. For all p > m = degree of the spline, the infinite norms 
of the spline dQIs Q* 2 are uniformly bounded by ^zj- This bound is 
independent of p and of the sequence of knots T. 

Proof: We have to find a good upper bound of 

v: = \K(-p)\ + \K(o)\ + \K(p)\ 
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where, expanding the determinants, we have 

Ki-p) = -^/Wi+p - Oi-vWi - Bi-p) 

\m = ^ + of ) /{Oi+p-OiWi-Oi-p) 
K(p) = -oi 2) /(9 i+P - e^Wi+P - Oi)- 

We recall that Oi = ^ Y1T=g h-r an d 

fl(2) _ 1 ( , \2 _ 

i ~ 2/ t\ / j \ L i—r L i—s) — n/ i\ 

m z (m — l) ^-^ m z (m — l) 

We first compute 

1 1 P 

Oi — Oi-p = — ^ (U-r — U-p-r) = — ^ ^ hi-k-r+1 = S 2 /m. 
reN m r£N m fc=0 

1 1 P 

#i+p — Oi = — (^i+p-r — *i-r) = — ^-i+fc-r+l = S^/m, 

reN m r6N m A;=0 

The proof being essentially the same for all p > m and for all i G K, we 
can restrict our study to the cases p = m and % = m — 1. In that case, we 
get 

m— 1 

£2 = mhi- m+ i + ^2 k(h + h) > S' 2 = h 1 + 2h 2 + . . . + (m — l)/i m _i, 
fc=i 

m — 1 

S3 = m^_ m+ i+y^ k(h+h) > S' 3 = (m-l)h 1 + (m-2)h 2 +. . .+2h m _ 2 +h m _\. 
fc=i 

Denoting, for 1 < /c < m — 1 

Sfe = ^i + • • • + hk, 

and s = s m -i, we get 

S 2 = si+s 2 +. • .+s m -i) S3 = s+(s-si)+(s-s 2 )+- • • (s-Sm-2) = ms-S 2 
whence 



S 2 ^3 > = ms(si + s 2 + • • • + s m -i) - (si + s 2 + . . . + s m _i) 2 . 
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Now, we come back to Si and we shall prove that Si < S' 2 S' 3 < S2S3. Si 
can be written under the form 

m — l m — l m — l m 

j=l i=l j = l i=j+l 

from which we deduce 

m — l m — l m 

Si = (m-l)Y,s 2 t -2j2 *j E *• 

i=l j = l i=j+l 

Moreover, for all 1 < i < m — 1, we have 

(m - l)s i = ms i - Si < msiS m -i - s { 
therefore, we obtain the result 

m — l m — l m — l 

Si<(m-l)J2 s 1 < ms m-i J2 s *~lL^ = S 2 S 3 ^ SA 

i=l i=l i=l 

Finally, for alH G TL 1 we have 

* 1 2 Si 2 ^ + 1 

v i = 1 + - 1 + T = 7' 

m — l m — l m — l 

whence HQ^IU < max ieZ i/f < ^±±. □ 

In the next section, we prove that the QIs Q* 2 are near-best in the sense 
of section 6. 



§8. A family of near-best spline discrete QIs 

For dQIs <5 P ,2 depending on p > m parameters, the coefficients (see 
proof of theorem 5) are given by 

K(-p) = -oi 2) /(e i+p - 9i- P Wi - 0i- p ) 

A*(0) = 1 + ei 2) /(9 i+p - 0i)((0< - (9i_ p ) 

a?(p) = -ol 2) /(e i+p - e t . P w l+p - 6i). 

Now, let us write the minimization problem of section 7 in Watson's form. 
Denote 

A; = (Xi(-p + 1), . . . , Ai(-l), Ai(l), . . . , lambdaiip - 1)) T e R 2p - 2 
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K = (AJ(-p), 0, . . . , A*(0), 0, . . . , \*(p)f G H 2 ^ +1 

Let G ]r( 2 p+ 1 ) x ( 2 p- 1 ) be the matrix with the following coefficients 
(notations of section 7) 

For r G K 1 : A^-p, r) = a r , A^O, r) = r , A^p, r) = ~7 r , 

For s G K 2 : Ai(-p, s) = -a s , Ai(0, s) = f3 8 , Ai(p, s) = 7s , 
For r G K 1 : A^r, r') = 0, r' ^ r, A^r, r) = -1, A^r, s) = 0, s G K 2 , 
For s G K 2 : A;(s, r) = 0, r G K±, A^s, s) = -1, A^s, s') = 0, s' ^ s. 
Then, using these notations, we can write 

ll-MIl = \\K ~ 

Theorem 6. Assume that the sequence of knots T satisfies, for all i E7L, 
the following properties 

+ 0i < 9i-p + &i+p < #i + 

then, for all i E 7L, X* is an optimal solution of the local minimization 
problem min || A^ || i . Thus, for all p > m, the spline dQIs Q* 2 are near-best 
and their infinite norms are uniformly bounded by ^-jy- This bound is 
independent of p and of the sequence of knots T. 

Proof: According to Watson's theorem, we must find a vector v* G R 2p+1 
satisfying 

IKIloo < 1, Ajv* = 0, v*(r) = sgn(X*(r)) for r = -p,0,p. 
Let us choose 

v*(-p) = -l,v*(0) = l,v*(p) = -l, 

V *(r) = -OL r + (3 r + 7 r, for T G K U 

v*(s) = -a 8 + (3 S +7 S , forseK 2 . 
Then it is easy to verify that the equations Af v* = are satisfied. More- 

—(2) 

over, the above expressions of A*(r) for r = — p, 0,p with 9\ ■ > imply 
that sgn(v*(r)) = sgn(X*(r)) for r = —p,0,p. It only remains to prove 
that, for (r, s) G K\ x K 2 

\v*{r)\ = | -a r + (3 r + 7r | < 1, \v*(s)\ = \a s +(3 s - ls \ < 1. 
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As (3 r = 1 — a r + 7 r for r E Ki and /3 S = 1 + a s — 7 S for s G K2, it is 
equivalent to prove 

< a r - 7 r < 1, < 7 S - a s < 1, /or (r, s) e K x x K 2 

We only detail the proof for r G Ki, that for s E K 2 being quite similar. 
Using the Vandermonde determinants, we get 

OL T — lr = Vi — r )(6 i+p — 9i-p)[(9i +p + Oi-p) — (0 r + 6i)], 

As 6i — 9 r > and 9i +p — di- v > 0, we shall have a r — 7 r > if and only 
if 

9 r + 9i < 9i +p + 9i- p 

for all r G -ftTi. However, since we have 9 r + 9i < 9i-\ + 9i, there only 
remains the unique condition 

9i-\ + 9i < 9i- p + 9i+ p . 

The other inequality a r — 7 r < 1 can be written 

{9% — 9 r )[(9i +p + 9i- p ) — (9 r + 9i)] < (6i — 9i-p)(9 i+p — 9i) 

Setting 8\ = 9 r — 9i- p , S 2 = 9i — 9 r , and £3 = 9i +p — 9i, the latter inequality 
can be written 

S 2 (S 3 - Si) < S 3 (S 2 + <5i), or 6^62 + S 3 ) > 

which is obviously satisfied. For s G K 2l the inequalities < 7 S — a s < 1 
are satisfied if and only if 

9i- p + 9i +p < 9i + 9i+±, 

whence the conditions on the sequence of knots. □ 

Remark. Theorem 6 imposes some conditions on the sequence of 
knots. For quadratic splines, we have studied arithmetic and geometric 
sequences: in both cases, the higher is p, the stronger are the conditions 
and for p — > +00, T is closer and closer to a uniform sequence. 
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9.1. Approximation of functions 

When a spline dQI Q is uniformly bounded independently of the partition, 
we can apply a classical result in approximation theory (see [5], Th 22, 
and [11], chapters 2 and 5): 

||Q/-/||oo<(l + ||Q||oo)doo(/,5) 

where S is the space of splines. In particular, when Q is exact on the 
space P m , then for / e C m+1 (I), one has 

||Q/-/||oo<C7r +1 ||/ m+1 |U 

for some constant C which does not depend on the given partition. There- 
fore spline dQIs give the best possible approximation order. More detailed 
results on error bounds are given in [1],[13] and [23]. 

9.2. Quadrature formulas 

Approximating Jj f by Jj Q^if), where Q^if) is the quadratic spline dQI 
of section 8, gives rise to an interesting quadrature formula 

f i QUf) = ^1^0 + ^2 ^(/) l Bt + fn+1 l Bn+1 

As it is well known, £ B = ^, J z B n+1 = ^ and J] Bi = 
for 1 < i < n. This formula is exact on P 2 , but in the case of a uniform 
partition, (see [21]), it is exact on P3 and provides an interesting comple- 
mentary formula to Simpson's rule in the sense that, in general, errors for 
both formulas have opposite signs. This will be detailed in another paper, 
together with applications to integral equations. 

9.3. Pseudo-spectral methods 

One can approximate the first derivatives of a given function / at the data 
sites 

O n = {do = 9i = ^(U-i + U), for l<i<n, 6 n+1 =t n }. 
by the derivatives of the quadratic spline dQI of section 8.1 

n 

Q* 2 f = f(t )B + + f(t n )B n+1 . 

i=i 

For interior points $i, 3 < i < n — 2, we obtain the general formula 
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which can also be written, by setting /_,■ = f(0j): 

(QZfYfti) = 7-{-<7i a i-lfi-2 + + {°i - cr- +1 )aj]/i_i 

lli 

[-a i c i _i + (a i -a- +1 )6 i + a i+ i]/ i + [(cr,-cr- +1 )c, + 6 i+ i]/ i+ i + a- +1 c i+ i/i +2 } 

For the first indices < i < 2, the coefficients are modified according to 
the convention ho = 0, which gives erg = 0, cr' Q = 1, ui = 1 and = 0. We 
thus obtain 

(QS/Wo) = |-{(ai - l)/o + 61/1 + ci/ 2 } 

(QWtfi) = ^-{(wi - l)/ + [<r 2 6i + a 2 a 2 ]/! + 
[a 2 ci + a 2 6 2 ]/ 2 + a 2 c 2 / 3 } 
(Q 2 /)'(^2) = ^{-o"2ai/o+[-^2&i + (c^2- si(/ma3)a 2 ]/i+[-o-2Ci+(o-2-<T3)6 2 +f3 a 3]/2 

^2 

+ [(a 2 - a 3 )c 2 + cr' 3 bs]f 3 + (73C3/4} 

In the same way, for the last indices n — 1 < z < n + 1, the coefficients 
are modified according to the convention ft n +i = and we obtain similar 
formulas for (QS/)'(0„_i), (QUY^n) and (Q^/)'(^ n+ i). In the case of a 
unifom partition, the formulas are given in [21]. Concerning error esti- 
mates, it is rather easy to verify that (Q^fYi^i) — f'i^i) = 0(h 2 ) where 
h = maxi<j< n /ij. A more detailed study will be done elsewhere. These 
results can be used in pseudo-spectral methods, as described for example 
in [12] and [29]. 
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